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THEORETICAL STUDY OF THE LATERAL FREQUENCY RESPONSE TO 
GUSTS OF A FIGHTER AIRPLANE, BOTH WITH CONTROLS 
FIXED AND WITH SEVERAL TYPES OF AUTOPILOTS 
By James J. Adams and Charles W. Mathews 


SUMMARY 


A theoretical approach has been used to determine the lateral fre- 
quency response of a fighter airplane to side gusts and to rolling gusts 
at a Mach number of 0.7 and an altitude of 30,000 feet. The frequency 
response and the power spectral density of the motion of the airplane in 
response to gusts were determined for the airplane with controls fixed 
and for the airplane in combination with three different basic types of 
attitude autopilots. 

The response to gust inputs for the airplane with controls fixed 
exhibited a large resonance associated with the Dutch roll mode of the 
airplane. When the airplane was combined with the various autopilots, 
the addition of yaw damping greatly reduced the resonance. The addition 
of autopilot components to supply heading stability and roll stability 
provided good regulation in that the yaw and roll responses to gusts of 
the airplane were greatly reduced. Autopilots that controlled side force 
to low values and that provided good course response to command signals 
resulted in large roll response to side gusts. 

A simplified transfer -function method of analysis for determining 
the response of the aiiplane -autopilot combinations was tried, and a 
comparison of the results of this method and the complete, three-degree- 
of -freedom method of analysis is given. It appears that the simplified 
method could be used for gust studies . 


INTRODUCTION 


One of the factors that affects the precision of such flight opera- 
tions as landing and air gunnery is the response of the ai r plane to gust 
disturbances. Therefore, it is highly desirable to study the motions of 
airplanes and airplane -autopilot combinations which result from gusts and 
to establish which configurations are best suited to min imi ze the motions . 
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The present paper is concerned with the theoretical lateral frequency 
response and the power spectral density of the motion in response to 
gusts of a fighter airplane both with controls fixed and with three 
different types of attitude autopilots - The purpose of the paper can 
be stated as: 

(1) The development of equations expressing airplane lateral fre- 
quency response to gusts . The method used follows that developed in 
reference 1 which describes the longitudinal response to gusts. 

(2) The examination of the lateral frequency response to gusts and 
the power spectral density of the lateral response to atmospheric tur- 
bulence in an attempt to evaluate the relative merit of the autopilots 
in controlling the airplane . 

( 3 ) The investigation of various simplified transfer-function 
methods of determining the response to gusts of the airplane -autopilot 
combination and the establishment of their merit. 

The results are presented as plots of the frequency response and 
power spectral density of the motion of the airplane, both with controls 
fixed and in combination with the autopilots, in response to side-gust 
and rolling-gust inputs. 


SYMBOLS 


n 


wing span, ft 

lift coefficient, — ^ - — 

|pSV2 


rolling -moment coefficient. 


yawing-moment coefficient. 


... LI 

|pSv2b 

N 

JpS^b 


Cy side-force coefficient, = — =L-^- 

JpSV 2 

D differential operator, d/d s 

acceleration due to gravity, ft/sec^ 


g 
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K autopilot gain 

Kg nondimensional radius of gyration about X-axis, kg-yh 

nondimensional radius of gyration about Z-axis, 

Kgg nondimensional product-of -inertia factor 

k X radius of gyration of airplane about X-axis, ft 

kg radius of gyration of airplane about Z-axis, ft 

L lift, lb 

L' rolling moment, ft -lb 

i x nondimensional tail length parallel to X-axis, based on wing 

span 

1 2 nondimensional distance (parallel to Z-axis) measured from X-axis 

to center of area of vertical tail 

m mass of airplane, slugs 

N yawing moment, ft-lb 

n frequency, cycles/sec 

p rolling angular velocity, radians/sec 

r yawing angular velocity, radians/sec 

S wing area, sq ft 

s distance measured in spans, tv/b v. 

t time, sec 

V velocity of airplane with respect to still air, ft/sec 

v velocity along Y-axis with respect to still air, ft/sec 

velocity of side- gust (positive gust produces sideslip 
to right), ft/sec 

Y side force, lb 
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a angle of attack, radians 

P total sideslip angle, p 0 + pg, radians 

P Q angle between X-axis and velocity vector V, radians 

Pg angle of sideslip due to side gust, radians 

S & total aileron angle, radians 

8 r rudder angle, radians 

K course, radians 

p relative-density factor, m/pSb 

p air density, slugs/cu ft 

cr sidewash angle, radians 

0 angle of roll, radians 

D0 nondimensional rolling gust, radians/span 

O 

heading or yaw angle, radians 

a* nondimensional circular frequency, 2jtnb/v, radians/span 

Subscripts : 

a, r denote aileron and rudder servos, respectively 

c command signal 

eg center of gravity 

t vertical tail 

wb wing-body combination 

e error signal 

\Jr, 0, p denote yaw, roll, and side -force information of feedback in 
autopilot loops, respectively 



MCA TN 3603 


5 


Dot over quantity indicates differentiation with respect to time. 

Stability derivatives are indicated by subscript notation; for 
example, 





dCn 

a 2* 

2V 

j *h_ 
a ^ 

2V 


The stability system of axes is shown in figure 1. 


DESCRIPTION OF AIRPLANE AND AUTOPILOTS 


All calculations made in this paper are for a present-day fighter, 
jet-propelled airplane with unswept wings. A drawing of the airplane 
is shown in figure 2. This particular airplane is used in this study 
because frequency responses obtained from flight tests are available to 
check the calculated frequency responses . The derivatives and mass 
characteristics of the airplane are listed in table I. The stability 
derivatives are obtained from unpublished data. 

The three types of autopilots which are studied in combination 
with the airplane can be described as follows : 

(1) An autopilot designated the type 1 autopilot combines yaw 
damping (rudder moves in proportion to a yawing -velocity signal), 
heading stability (rudder moves in proportion to a heading error signal), 
and roll stability (ailerons move in proportion to a roll-angle error 
signal). This configuration was chosen for study because it represents, 
in an idealized form, a type of regulator autopilot that is being used 
at the present time. 

(2) An autopilot designated the type 2 autopilot combines yaw 
d amp ing, roll stability, a signal from the side fofce to the rudder 
(rudder moves in proportion to the side force to reduce the side force), 
and a signal, from the heading to the ailerons (ailerons move in proportion 
to the heading error to correct the heading error by rolling the airplane). 
This configuration was chosen for study because it represents a type of 
fighter autopilot being used at the present time . 
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(3) An autopilot designated the type 3 autopilot combines yaw 
damping, hea ding stability, roll stability, and a signal from the side 
force to the ailerons (ailerons move in proportion to the side force to 
tilt the lift vector in the direction of the side force. This configu- . 
ration was chosen for study because it has good response in both heading 
and course when a command is applied to the heading loop. 

The autopilot loops which the different autopilots have in common 
are given the same gain in each case. These loops, the yaw-damping 
loop and the roll-stability loop, are given the gain recommended by the 
manufacturer of a type 2 autopilot for use with an airplane of the type 
studied in this investigation. Since there are no other autopilot loops 
placed inside the yaw- damping and the roll-stability loops in any of the 
autopilots studied that may change their optimum gain setting, it is 
reasonable to assume that they will have the same gain for each of the 
autopilots . The heading-stability-loop gain of the type 1 and type 3 
autopilots is adjusted to give a closed-loop heading response to rudder 
deflection having a damping ratio of 0.3 (based on the assumption of a 
single degree of freedom in yaw with the yaw-damping gain mentioned 
previously) . 

One calculation with the type 2 autopilot is made with the gain 
relating the heading error signal to the aileron deflection and the gain 
relating the side-force signal to the rudder deflection set at the values 
recommended by the manufacturer of a type 2 autopilot. A second calcula- 
tion is made in which the heading-error-signal gain is increased so that 
the steady-state roll angle resulting from a heading error for the type 2 
autopilot will equal the steady-state roll angle resulting from a heading 
error for the type 3 autopilot. This gain adjustment affords a direct 
comparison between these two autopilots and, in addition, provides a 
rapid course response in each case. In the type 3 autopilot the gain 
between side force and aileron deflection is chosen to result in a roll 
angle of 1 radian for a side-force signal of 0.05g acceleration. All 
autopilot gains are listed in table U. 


DESCRIPTION OF GUSTS 


Two types of gusts are studied - side gusts, which add an increment 
to the angles of sideslip, and rolling gusts, which have the effect of 
adding an increment to the rolling velocity of the airplane with respect 
to the air mass. These gusts are assumed to vary sinusoidally. The use 
of sinusoidal disturbances appears to be appropriate for the study of the 
effects of turbulence*, since flight data obtained in rough air show that 
the disturbances are of an irregular oscillatory nature and since such 
irregular disturbances can be reduced by generalized harmonic analysis to 
sinusoidal components of various frequencies and amplitudes . In this 
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paper the disturbances are restricted to side gusts and rolling gusts 
inasmuch as reference 2 has shown that the lateral effect of turbulence 
can he well approximated by gusts of these two forms . These gust forms 
are convenient to use since they are expressed by analytical terms that 
are similar to the terms for the variables 3 and D0 ordinarily used 
in equations of motion. Furthermore ; references 2 and 3 show that the 
relative magnitude of the power spectral density of the side gusts 
diminishes as the square of the frequency and that the rolling gusts have 
a power spectral density that is constant throughout the frequency range 
considered in this paper. 

The frequency range covered in this paper varies from 12 radians 
per second down to 0.2 radian per second. This range extends to well 
above the primary resonant frequency of the systems studied. The upper 
limit represents a gust wave length that is still large with respect to 
the length and span of the airplane. The lower limit of 0.2 radian per 
second was chosen because it was believed that information concerning 
gust frequency response is of little importance at lower frequencies 
inasmuch as the pilot or autopilot can easily control the low-frequency 
lateral motions . However , the lower limit is too high to show completely 
the effects of the spiral mode of the airplane with controls fixed. 


METHOD OF ANALYSIS 

Equations of Lateral Motion Including Gust Inputs 


The equations of lateral motion used in this investigation are 
.similar to those presented in reference 4. The usual procedure is varied 
by expressing separately the force and moment contributions of the wing- 
fuselage combination and the tail. 

The equations written in nondimensional form are 

A 

* w ) • * C Yp t Pt + + |plp + i^L + C Y 8 ^ & r + [ 


2P - KxzD 2 ^) 


Svb PcS 




D0, 


eg 


+ C n + ic_ B0+ + cu a. + C n 


( 1 ) 


2^(k x 2 D^ - KxzD 2 *) = C X + If 


Pwi 
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In these equations the coefficients of the tail forces and moments 
/Cy , C n , etc A are based on wing area and wing span. 

\ 0t Pt ) 

In order to obtain the response to gusts from equations (l) the 
values of p cg , $t> 0 cg> ^cg mua- t be expressed in terms of 

the gust inputs and the airplane motion. Derivation of the expression 
for these angles requires consideration of the assumed form of gust 
disturbance and the sidewash effects at the tail. The gust disturbance 
is divided into two separate modes - side gusts and rolling gusts. 

Side gusts are assumed to consist of lateral sinusoidal, velocity varia- 
tions within the air mass at right angles to the plane of symmetry of 
the airplane. Rolling gusts are assumed to consist of rotating sinus- 
oidal velocity variations within the air mass and are assumed to rotate 
about the path of the airplane. At a given point in the air mass, the 
gust velocity is assumed not to change with time, and the wave length of 
the gust is assumed to be long compared with the length and span of the 
airplane . 

The equations are solved for v/V, and f in terms of the gust 
velocities. A side gust with velocity v g will change the angle of side- 
slip of the airplane an amount equal to v g ^Vj that is, p g = v g fv . In 

addition, the airplane will have an angle of sideslip relative to still 
air equal to v/V; g 0 = v/v. Then the total angle of sideslip at the 

center of gravity is 


Peg Po + Pg 

likewise, the total rolling velocity of the airplane may be considered 
as consisting of the rolling velocity of the rolling gust and the rolling 
velocity of the airplane in still air: 


Also, 


^cg = ^ 



= V 


The angle of 
angle of Bide slip 
velocity, and the 


sideslip at the tail is assumed to be a function of the 
of the airplane, the sideslip at the tail due to yawing 
sidewash angle resulting from wing lift which existed 



NACA TO 3603 


9 


when the wing was at a position now occupied by the tail. There is, 
therefore, a lag that occurs during the time required for a given lift 
disturbance at the wing to move back and cause some associated sidewash 
at the tail. This lag is proportional to the time required for the 
airplane to travel 1 tail length. Sidewash can be considered to arise 
from two sources: sidewash due to sideslip, both p Q and p g ; and 

sidewash due to rolling velocity, both D0 and D0 g . Estimations of 

sidewash magnitude are given in references 5 and 6. The effect of p g 

at the tail is assumed to lag behind the effect of p at the wing. 

S 

This lag describes the fact that a given gust disturbance is encountered 
at the center of gravity prior to being encountered at the tail . The 
expression for angle of sideslip of the tail expressed in nondimensional 
notation is 



The rolling velocity of the tail is assumed to equal the rolling velocity 
of the airplane and the rolling velocity of the gust, with lag to repre- 
sent the effect mentioned previously: 



D0 


KL-z, 


(3) 


The constant time lag -l x is expressed as 1 - I x D, which are the 

first two terms of a power-series expansion for a constant t ime lag. 

Using the expression 1 - Z x D to denote lag in sidewash has been shown 

to agree very well with experimental measurements of lag in sidewash. 
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Including these lag effects is an improvement over previous calculations 
which have neglected these effects . For a fighter airplane in the flight 
conditions used in this paper, the lag effects are negligible . The equa- 
tions are set up in complete form for possible use where the effects may 
be more important. The lag in sidewash terms may be important for air- 
planes with sweptback wings or at flight conditions requiring high lift 
coefficients . The lag in gust effect at the tail is important for large 
airplanes, airplanes at low velocity, or airplanes with low-density 
ratios. These terms can be included with very little added complication. 

The expression for angle of sideslip of the tail and the rolling 
velocity of the tail can now be written as follows: 

H - So - + Pgt 1 - D ) - Polfl 1 - D ) - 

D ) - D ) ■ D ) {k) 

and 

D0 t = D0 + D0 g (l - l x D) (5) 

Substituting the expressions for p cg , D0 cg , and into 

the equations of motion (eqs. (l)) gives 

+ 2^ - %Po - - 

% p e + l(% + c 0 1 *8 + + ' ^t}^ 8 

- w° - " 2 V* " is 130 ■ ^ ‘ Cn Ba Sa “ 

(6) 

S P e + KSp + t Vr) 1W B + iV^ + 

SMKX 2 ^ ~ 2UKXZ D ^ - C l p Po * _ g 0 ^ 1 ^ " ' Cl 5 r Sr “ C *s a 8a = 


% p e + l^inp + c irJ Dp s + + ‘ Ci dp 1 ) d2?, 8 
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■where 


C Yq + f 1 - If) 

Pvt> Pt \ d P/ 

% ■ %*, * vt 1 ■ t) 

% = C ^ b + sj 1 - if) 


%s ■ 

° l Dp “ 2I *if ° i 3 t 


% ■ 

° ir = Cl r vb - 2 ^^ t 
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n P = 


C, 2 [ 




C, - 2l x &L-C, 
*D^jS SD0 l 3 t 



2 : x x zCy 

Pt 


Sp t " 2VzGn 3 t 



2 Wi. 
Pt 



- -^Pt 


Note that 



2Z 




2 ^ c n, 


Pt 



2^1 


Pt 


The terms involving Df3g and D^0g on the right-hand side of equa- 
tions (6) result from the lag in sidewash and from the lag between the 
effects of a gust on the wing and on the tail as the airplane passes a 
given point in the air mass . 
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Calculation of Response to Sinusoidal Gust Disturbances 

Airplane •with, controls fixed .- To determine the steady-state response 
of the airplane to a sinusoidal gust disturbance when the controls are 
fixed, terms involving 5^. and S a are set equal to zero. The equations 

of motion can then be solved, simultaneously (for example , by the method 
of determinants) to obtain the variables p Q , ty, and <f) in terms of 

p and. D0~. The steady-state response to a sinusoidal input can then be 

CD O 

found by substituting D = its into the resulting expressions , which are 
known as the transfer functions . 

Airplane -autopilot combinations .- A second aspect of the airplane 
gust response studied, was the motions of the airplane encountered when 
autopilot loops were incorporated in the airplane to regulate some of the 
variables. Of the many methods of analysis available for dete rmin i n g the 
response of an airplane -autopilot combination, two methods were used: 
an approach using the equations of motion and an approach using component 
transfer functions . For each calculation the autopilot servomechanisms 
were assumed to have perfect response (no lag and constant amplitude) over 
the range of frequencies of interest in the calculations of the airplane 
response. Therefore, the transfer functions of the autopilot servo were 
expressed as gains. 

In the equation-of -motion method of calculation, the effects of the 
autopilot loops were included as equivalent stability derivatives. This 
substitution was possible because the assumption of a perfect servo makes 
the forces and moments applied by the autopilot on the airplane explicit 
functions of the motion of the airplane. A three -degree -of -freedom 
response for the airplane including these artificial stability deriva- 
tives was obtained. The analysis was not si gnif icantly more complicated 
than that required to obtain the response of the airplane with controls 
fixed. calculation was repeated for the airplane in combination with 

each of the three autopilots described previously. The equivalent sta- 
bility derivatives are listed in table II beside the autopilot gains from 
which they were derived. 

The transfer-function method of analysis also was used in this 
investigation to present examples of the application of this type of cal- 
culation in the determ i nation of the response to gusts. This method is 
very useful to engineers in establishing optimum autopilot gains, since 
it affords a direct study of the effects of gain change on the system. 

The method is also used for cases in which the operation of the autopilot 
cannot be expressed analytically. Many times in such cases the actual 
autopilot equipment is used in conjunction with the analog-computing 
equipment . 
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The transfer-function method is carried out by combining the transfer 
functions of the airplane with the transfer functions of the autopilot in 
a manner that is discussed in reference 7» A block diagram illustrating 
this method is shown in figure 3- In general, the transfer function of an 
actual autopilot servo would be frequency- variant and would be expressed 
by a function of second order, or higher. However, as was stated before, 
the autopilot servos are assumed to be perfect since study of the effect 
of variations of "servo performance was considered beyond the scope of the 
present investigation. 

Even with this simplified assumption for the autopilot servos, the 
combination of a three -degree -of -freedom airplane with an autopilot 
according to the scheme shown in figure 3 will result in a system trans- 
fer function having polynomials of approximately the fourteenth order. 

If servo dynamics were included, the order of the transfer function 
would be even higher. 

Since it is very difficult to study parameter change with tr ans fer 
functions of this complexity, various simplifications were tried in an 
attempt to find the simpliest form of system transfer function that would 
give an answer consistent with the results obtained from the analysis of 
the equations of motion previously discussed. Details of these simplifi- 
cations are given in the appendix. 

A comparison of the equation-of -motion method of analysis with the 
transfer-function method is interesting. The equation-of -motion method 
using the artifical stability derivatives and the transfer-function 
method describing a three -degree -of -freedom airplane and using the sim- 
plified assumptions for the autopilot servos both express the same system; 
the only difference is that the equation-of -motion method results in a 
system transfer function of the fifth order and the transfer -function 
method results in a system transfer function of the fourteenth order. 

These facts lead to the conclusion that there are factors in the numera- 
tor and denominator of the system transfer function obtained by the 
transfer-function method that could be extracted so as to reduce the 
order. However, extraction of these factors is too difficult to be 'useful 
in simplifying the method. 


Calculation of Power Spectral Density of Airplane Motion 

The power spectral density of the side gusts (3g and the rolling 
gust D0g were derived from reference 2 and are plotted in figure 4. 

The square of the frequency response to p g and D0 g of the motion of 

the airplane as obtained by the equation-of -motion method was multiplied 
by the corresponding power spectral density of the gust to determine the 
power spectral density of the motion of the airplane. These curves have 
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been stopped at the low -frequency limit at which the variation of the 
power spectral density for the side gusts appears to break. Recent 
information indicates that this falling off occurs at gust wave lengths 
of approximately 6,000 feet. 


Calculation of Frequency Response to Command Inputs 

General . - As a point of additional interest, the heading response 
to a heading command and the course response to a heading command for the 
airplane -autopilot combinations are also- calculated. The command response 
considered in this paper is one that results from applying the command to 
the particular autopilot channel that is vised to regulate the variable 
being commanded. In the case of pilot-applied heading or course commands 
for the type 1 autopilot, this method is not usually used because such 
commands are applied to the aileron loop while the heading regulation 
supplied by the yaw loop is rendered inoperative during the command oper- 
ation. When commands result from tying the airplane to an external ref- 
erence, such as radar tracking of a target or beam following, it is not 
possible to utilize separate loops for command and regulation. Systems 
in which the same loop is used for co mman d and regulation are considered 
herein. 

Type 1 autopilot .- The heading response under the foregoing condi- 
tions for the type 1 autopilot may be determined as the response to a 
rudder deflection, and the course response may be determined as the 
response of the path of the airplane to a rudder deflection. The heading 
response ijr/ijr c is determined by using a rudder input in the equations of 

motion, similar to the procedure used to determine heading response to 
side gusts t/Pg* The expression for the course response to a heading 

command is developed as follows : 


« = t -k P 0 
k _ t ^0 

*c 


= JL + JL 

t c V t c 
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where p 0 /t is the ratio of the quantities p 0 /&r and \|r/6 r , which were 
obtained from the equations of motion. 

Type 2 autopilot .- In determining the response to a heading command 
for the type 2 autopilot, the heading signals go to the aileron loop. 

With the lateral force held to small values by the side-force component 
of the autopilot, the resulting heading and course c hang es are the same. 
These responses are obtained by assuming that the horizontal sideways 
force, horizontal acceleration, and consequently the rate of change of 
heading and course are proportional to the roll angle 0 . The horizontal 
sideways acceleration is integrated to obtain the change in heading or 
course. The development of the equations of these relationships is as 
follows . The open-loop response of 0 to the error signal is 


H e 


0 

1 + K0 f-Ka 

Uq 


The relationship between t and 0 is 


ir = 

V 2 D 


Therefore, the closed-loop response to a heading co mman d is 


jr _ k 

H c *c 


¥ o. 


' a gp 

1 + K ?5 |-K a V 2 D 

£L 


1 + 


0 


1 + 


0 


gb 

V%) 


(7) 
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Type 3 autopilot .- The heading response \|r/i|r c for the airplane in 

combination with the type 3 autopilot is obtained in the same manner as 
for the type 1 autopilot. In order to determine the course response to 
a heading command, the open-loop response of 0 to the error signal 

is obtained by including a signal in the roll equation that results from 
the coupling of side force to aileron deflection of the autopilot. Ibis 
signal, which is similar to an equivalent stability derivative, is given 
in table II. The resulting value of 0/^e is treated in the same way 

as for the type 2 autopilot case . The horizontal sideways acceleration 
was assumed to be proportional to 0, and this acceleration was integrated 
to give the change in course. This open-loop response was closed with a 
unity feedback to give values of tc/i|r c . 


RESULTS 

Frequency Response of Airplane With Controls Fixed 


The response to side gusts and to rolling gusts of the airplane with 
controls fixed is shown in figure 5 . These results show that, at low 
frequencies, the yaw response to side gusts ty/Pg is approaching a value 

of 1.0; this trend indicates that the airplane yaws into the gust. The 
same result can be stated another way by saying that p is approaching 
zero at low frequencies. The roll response to side gusts 0/p g is pri- 

BKirlly the result of the C lf> couplihg. . A large resonance occurs at the 

frequency of the Dutch roll oscillation of the airplane. The same 
resonance appears in the response to rolling gusts . These resonance 
peaks are due to the low damping of the Dutch roll mode, a characteristic 
that is fairly typical of present-day fighter airplanes . Moreover, the 
roll response to side gusts is larger than the roll response to rolling 
gusts at the peak frequency whereas the comparison is reversed at higher 
and lower frequencies because the Dutch roll mode is excited more by side 
gusts than by rolling gusts . 


Frequency Response of Airplane -Autopilot Combinations 

Type 1 autopilot .- The frequency response of the airplane in combi- 
nation with the type 1 autopilot is shown in figure 6. When the airplane 
was combined with the type 1 autopilot, the resonant peaks that were 
present in the responses of the airplane with controls fixed were greatly 
reduced, and the amplitude ratio of the t/Pg responses was reduced by a 

factor of approximately 10 in the low frequency range. The additional 
heading stability provided by the autopilot approx i mately doubled the peak 
frequency of the i{r/pg response over that of the airplane with controls 
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fixed. Since the airplane was restrained from heading into the gust, 

P was increased beyond the values for the condition of fixed controls, 
and, in spite of the fact that the autopilot stabilized the airplane in 
roll, an increase in the 0/Pg response resulted at low frequencies. The 

roll response could be reduced by increasing the roll loop gain. The 
roll gain used in this calculation is relatively small, and, with some 
types of autopilot servos in present-day use, a higher gain could be 
used without encounter in g objectionable characteristics. The roll sta- 
bility added by the autopilot produced a new resonant peak for the roll 
response at a frequency which is not very different from that of the peak 
of the airplane with controls fixed due to the Dutch roll oscillation. 

The roll response to rolling gusts 0/ D0g was greatly reduced at low 

frequencies by the autopilot. 

The results of calculations of a simplified block di agram also are 
plotted in figure 6 for comparison with the results of the three -degree - 
of -freedom equation of motion. The amplitude ratio and the phase angle 
of the yaw response to rolling gusts \|r/D0 g at low frequencies for the 

simplified calculation do not agree with those for the three-degree-of- 
freedom calculation. The predicted amplitude in both cases is small, 
and this fact is probably the reason for the poor agreement. Although 
in some applications of frequency-response techniques - for example, in 
the study of the open-loop response of a mechanism - such a phase -angle 
error would be unacceptable, the phase relationships are not believed 
to be particularly significant when gust inputs are considered. 

In addition, the control-deflection frequency response as calculated 
by the transfer -function method is also plotted in figure 6 . The figure 
shows that there are only small variations in amplitude ratio and phase 
angle with frequency up to the peak frequency. At higher frequencies the 
inertia effects attenuate the airplane response, and consequently the 
control deflections are reduced also. 

Type 2 autopilot .- When the airplane is combined with the type 2 
autopilot with the lower value = 23 volts/radian) for the autopilot 
coupling between i|r and Sg,, the response to a side gust, as shown in 
figure 7 , indicates that the resonant peaks have been greatly reduced. 

The low-frequency values of \|r/p g are the same as for the airplane with 

controls fixed; the similarity is to be expected since the side-force 
component of the autopilot aids the airplane in turning into the gust. 

The values of 0/ pg are increased at the low frequencies by the auto- 
pilot coupli n g between i|r and 0. The plot of the response to rolling 
gusts (fig. 7 (b)) shows that the incorporation of roll stability has 
reduced the 0/D0 g response at low frequencies. The i{r/D0g response 
is small, although it is somewhat greater than that obtained with the 
type 1 autopilot, particularly at lew frequencies . Once again the 
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results of the simplified block-diagram method are plotted with the 
results of the three -degree -of -freedom equation of motion for comparison. 
It can be seen that some discrepancies occur in the two calculations for 
\)r/D0g. The control-deflection curves show a rapid increase in rudder 

amplitude ratio with increase in gust frequency. The aileron amplitude 
ratio for side gusts follows the variation of ijr/Pg and has the same 

phase angle. Both the variations of 0/D0 g and S^I^g as calculated 
by the simplified method are exactly the same as for the type 1 autopilot. 


When the autopilot coupling between ijr and 8 a is increased 
= llA volts/radian), the 0/p g response is also increased. (See 


fig. 8.) This coupling, discussed subsequently, wan increased to pro- 
vide better heading response to a heading command, and this heading- 
response improvement is obtained at the expense of a higher amplitude 
ratio for 0/p g . The effect of the large roll angles that occur at low 

frequencies with the type 2 autopilot causes a reduction in the ijr/p g 
response. The simplified block-diagram method of calculating ilr/p g , 

0/D0g, and i}r/D0g is the same as for the case with Ify = 23 volts/radian. 
The effect of the large roll response to side gusts, which was ignored in 
the simplified calculations, results in discrepancies between the simpli- 
fied solutions and the equation-of -motion solutions at low frequencies. 


Type 3 autopilot .- The results of the airplane in combination with 
the type 3 autopilot are shown in figure 9* The coupling between side 
force and aileron deflection, which was included to provide good course 
response to a heading command, caused a large roll response to side gusts 
to occur. The f/Pg response is the same as for the airplane in combi- 
nation with the type 1 autopilot with the exception that the large roll 
angle that occurs at low frequencies with this configuration causes a 
reduction in the t/Pg response. The response to rolling gusts is, in 

general, of the same order of magnitude as for the other combinations . 


The simplified block-diagram solutions for the airplane in combina- 
tion with the type 3 autopilot are in poor agreement with the equation-of- 
motion solutions. The simplified solution for ty/Pg is in error because 
it- did not account for the effects of the large roll angles that occur. 

The solution for 0/p g is in error because the simplified expression 

for ty/Pg is included in the expression for 0/p g ; and, because of the 
large coupling between side force and roll angle, the solution of 0/Pg 
is very sensitive to any errors in .t/pg- Thus, even the minor errors 
in the simplified ijr/p g response at high frequencies are reflected as 
large errors in the 0/Pg response. A modified solution for 0/ p g , in 
which a three-degree-of -freedom \jr/p g response, where the yaw loops of 
the autopilot were included as equivalent stability derivatives, did show 
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good agreement with the equation- of -motion solution. (See fig. 9-) The 
modified solution also demonstrated the reason for the unusual shape of 
the frequency-response curve for 0/Pg. The increase in 0/p g with 
frequency at low frequencies is caused by the increase in side force due 
to sideslip and the coupling between side force and roll angle. Above a 
frequency of 0.1 radian/span, the side force due to rudder deflection 
increases rapidly and with an opposite sign from the side force due to 
sideslip , until at co = 0.22 radian/span the two effects approximately 
cancel. At frequencies higher than = 0.22 radian/span, the side 
force due to rudder deflection and the side force due to sideslip become 
in phase; the result is an increase in side force and consequently in 
0/Pg. This increase continues until the inertia effects attenuate the 
response . 

The simplified solutions for \(r/D0g and 0/ D0g are identical to 

the solution for the type 1 autopilot, but are in only fair agreement 
with the equation-of -motion solutions . A modified solution for 0/D0g 

similar to that described for 0/3g, also shown in figure 9> gave a 
response that was in good agreement with the equation-of-motion solution. 


Power Spectral Density of Airplane Motion 

The power spectral density of the motion of the airplane in response 
to gusts is shown in figure 10. These plots emphasize the greater rela- 
tive importance of the low-frequency side -gust components over the high- 
frequency components and the greater relative importance of the. side gust 
over the entire frequency range compared with the rolling gusts in pro- 
ducing motions of the airplane. The area under these curves over any 
particular frequency is proportional to the mean square of the part of 
the airplane motion contributed by those frequencies. Aside from these 
general conclusions, the power-spectral-density plots also show the vari- 
ations between configurations that were noted previously in the discussion 
on the frequency response . These similarities existed in the cases con- 
sidered in this paper ina smuch as all comparisons were made for restricted 
frequency ranges. 


Summary of Results on the Response to Gusts 

When the airplane was combined with the various autopilots, in all 
cases the Dutch roll resonance was greatly reduced. The heading stability 
provided by the type 1 autopilot maintained good regulation in yaw, but 
the increase in sideslip at low frequencies resulted in a moderate increase 
in roll response to side gusts . This roll response could have been 
reduced if greater roll stability had been provided by. the autopilot. The 
response to rolling gusts was well regulated. When the airplane was com- 
bined with the type 2 autopilot, the airplane yawed into the gust to 
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reduce sideslip, and the coupling between heading error and aileron deflec 
tion caused large roll response to side gusts. The response to rolling 
gusts was approximately the same as for the type 1 autopilot. When the 
airplane was combined with the type 3 autopilot, good regulation in yaw 
was maintained. Ths coupling .between side force and aileron deflection 
caused large roll response to side gusts . The response to rolling gusts 
was well regulated. 


Response, to Heading Commands for the 
Airplane-Autopilot Combinations 

In order to complete more nearly the Information necessary for 
evaluation of the different autopilots discussed in this paper, the 
response to hea ding commands for each of the airplane -autopilot combina- 
tions is also given. The heading response to a heading command for all 
cases is plotted in figure 11, and the course response is plotted in 
figure 12. 

The excellent heading response to a heading command for the airplane 
in combination with the type 1 and type 3 autopilots is illustrated in 
figure 11. Each maintains an amplitude ratio of approximately 1 out to 
frequencies well beyond the Dutch roll frequency of the basic airplane. 
The heading response of the airplane with the type 2 autopilot is the 
poorest. However, increasing the coupling of | to ^ improves the 
heading response. The larger value of the coupling of \]r to 0 used in 
the calculations was chosen to give the same static ratio of roll angle 
to heading command as the type 3 autopilot. Therefore, the frequency 
response for these two cases can be compared directly. 

A comparison of the course response to heading comnnands for the 
three autopilots (fig. 12) indicates faster course response for the 
type 2 and type 3 autopilots than for the type 1 autopilot. The reason 
for this difference is that for the type 1 autopilot only side force is 
available for producing horizontal -path changes, since, as was explained 
before, in the present study the command was restricted to the yaw loop. 
For the type 2 and type 3 autopilots the lift force is tilted to provide 
a horizontal component of acceleration. The type 3 autopilot combines 
both the good heading response of the type 1 autopilot, by including the 
same heading-stability loop as the type 1 autopilot, and an excellent 
course response, by adding the coupling between side force and roll angle 
which tilts the lift vector. 

In many flight operations, it is desirable to have rapid course 
response to automatic lateral commands, such as for rocket-armed fighter 
operations or for beam riding. It is also desirable to regulate sideslip 
during pilot -commanded lateral maneuvers . However, the methods incor- 
porated to provide good course response and to regulate sideslip in the 
autopilot configurations considered in this paper result in large roll 
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response to side gusts. It might be possible to provide a dual mode of 
operation that would provide sideslip regulation at low frequencies and 
heading regulation at high frequencies. Since i]r = p 0 at high frequen- 
cies, this arrangement would provide sideslip regulation thro ugho ut t he 
frequency range insofar as command maneuvers are concerned and would 
avoid the large roll motions that occur because of high-frequency side 
gusts . 


CONCLUDING REMARKS 


In developing the equations of motion as used in this paper to 
express motion due to gusts, terms involving sidewash at the tail and 
expressing lag effect of the gust at the tail were used. However, for 
the airplane configuration vised in these calculations, these terms had 
little or no effect and could have been neglected. For airplanes with 
sweptback wings or for high values of lift coefficient, these terms may 
become significant. The lag term which represents the lag in gust effect 
at the tail may be important for large airplanes, airplanes flying at 
low speed, or aiplanes with low-density ratios. 

The response of the airplane with controls fixed exhibited a large 
resonance associated with the Dutch roll mode. This resonance response 
was particularly noticeable in the roll response to side gusts. This 
feature of a poorly damped Dutch roll mode that is excited by side gust 
is typical of present-day fighter airplanes. At low frequencies the 
airplane yawed into the gusts and reduced sideslip due to gusts. 

For cases in which various attitude autopilots were combined with 
the airplane, the following changes to the response of the airplane with 
controls fixed are noted: 

(1) When the airplane was combined with an autopilot that supplied 
yaw damping, heading stability, and roll stability, the increase in side- 
slip at low frequencies resulted in a moderate increase in roll response 
to side gusts. The heading stability maintained good regulation in yaw. 

( 2 ) When the autopilot supplied yaw damping, roll stability, a 
coupling between the heading error and the aileron deflection, and a 
coupling between side force and rudder deflection, the airplane .yawed 
into the gust in order to regulate side forc'e and, consequently, side- 
slip. The coupling between heading error and roll angle caused large 
roll response to side gusts . 

(3) When the autopilot supplied yaw damping, heading stability, roll 
stability, and a coupling between side force and aileron- deflection, the 
heading stability again maintained good regulation in yaw. The coup ling 
between side force and aileron deflection caused a large roll response 

to side gusts . 
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(4) In all three cases the yaw damping eliminated the resonant 
Dutch roll response to gusts. Also, in all three cases the roll response 
to rolling gusts was approximately the same, with the roll stability added 
by the autopilot maintaining good regulation in roll. 

In the last two cases, the type 2 autopilot with a heading-stability 
gain of 114 volts/radian and the type 3 autopilot, the features that pro- 
duced good course response to heading commands caused large roll response 
to side gusts to occur; this response was in contrast to the better roll 
response of the airplane in combination with the type 1 autopilot. 

The simplified method of calculation predicted the important motions 
of the airplane with good accuracy. Some discrepancies were noted in the 
predicted amplitude of the small motions . However, since the trends of 
the changes in the motion of the airplane with changes in autopilot con- 
figuration were accurately predicted and since the peak frequencies were 
accurately predicted, it appears that the simplified method could be 
used for gust studies. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 28 , 1955* 
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APPENDIX 

CALCULATION OF FREQUENCY RESPONSE OF AIRPLANE -AUTOPILOT 
COMBINATIONS BY TRANSFER-FUNCTION METHOD 


The details of the simplifications and the steps in calculating the 
frequency response of the airplane -autopilot combinations by the transfer- 
function method of ana lysis are described. The simplifications tha t were 
examined generally consisted in attempt ing to decrease the n umb er of 
'terms in equations (6) by reducing the degrees of freedom involved in the 
calculations . Motions other than the particular one being studied were 
ignored either because they were small or because they were controlled by 
some component of the autopilot and, therefore, did not influence the 
motion of the airplane in the mode being studied. The assumption that a 
particular motion can be ignored because it is controlled by some compo- 
nent of the autopilot is subject to the performance of that component. 

The loop gain must be large enough to regulate the variable to a point 
where it can be ignored, and the dynamics of the loop must be good eno ugh 
to provide this regulation throughout the frequency range being studied. 


Type 1 Autopilot 

For the airplane in combination with the type 1 autopilot, the 
responses in i|r, 0, S r , and 5 a to side gusts and rolling gusts were 

calculated. A block diagram of the autopilot setup is shown in figure 13. 
Since the airplane was controlled in roll, freedom in roll is ignored 
when calculating ljr/p g . Also, since the side force of the airplane was 
small in comparison with the yawing moments on the airplane, freedom in 
lateral displacement is ignored, and ijr is assumed to equal -p Q . On 
this basis, \jr/ p g and ty/Sr for the airplane with no autopilot yaw 

components are obtained by considering an airplane with a single degree 
of freedom in yaw. This airplane was combined with the autopilot compo- 
nents as shown in figure 13 to obtain the response of the controlled 
airplane . 

In order to obtain 0/p g for the controlled airplane, the strong 
coupling between 0 and p due to dihedral effect must be considered. 

In this case, since i}r is controlled, freedoms in yaw and lateral dis- 
placement are ignored, and both ijr and p Q are assumed to be eq ual to 
zero. Then, from the simplified roll equation (obtained from eqs. (6)) 
an expression for 0/p of the airplane with controls fixed is determined. 
If p is assumed equal to p g , the quantity 0/p g for the airplane 
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without roll control can be determined. This response of the airplane 
to Pg is added to the single -degree -of -freedom response of the airplane 
to 8 a in a mann er indicated by the block diagram of figure 13 to give 
the controlled-airplane response to Pg. The assumption that p = pg 
will be incre asing ly in error as the frequency is reduced because the 
airplane will drift with the gust at very low frequencies. However, at 
to = 0.01 radians/span this assumption is still accurate enough for 
practical purposes. 

In the dete rmina tion of the roll response to rolling gusts 0/D0 g , 
the freedoms in yaw and lateral displacement are ignored for the reasons 
given previously, and the response is obtained by considering an airplane 
with a single degree of freedom in roll in combination with the autopilot 
loops as shown in figure 13 (b). 

- In order to obtain the yaw response i]r to rolling gusts D0g, the 

couplings between 0 and i|r, 0 and p, and p and i|r must be con- 
sidered. Therefore, an expression for \{r/0 containing these couplings 
is determined from the simplified yaw and side-force equations (obtained 
from eqs. (6)). This expression is used to relate the homogeneous 
response of \[r to 0. This response in ijr is added to a single -degree - 
of -freedom yaw response to rudder deflection, in which the rudder deflec- 
tion is controlled by the yaw componenets of the autopilot, as is shown 
in figure 13 . The t/D0g response is then- obtained by multiplying the 

ty/0 response of the controlled airplane by the previously determined 
0/D0 g response . 


Type 2 Autopilot 

The various responses of the airplane in combination with the 
type 2 autopilot are calculated in the same basic manner as for the type 1 
autopilot combination except for the following changes in detail. Since 
the side force of the airplane is small and since the airplane is stabi- 
lized in roll by the autopilot, the roll and lateral displacements are 
ignored in representing the airplane in yaw with controls fixed. Again 
ijr is assumed to equal -p Q . -^ n or< 3- er to establish the effect of the 
autopilot side-force loop, the side force, even though it is small, must 
be calculated. From the equations of motion it can be seen that the 
nondimens iopal side force is 


o 


C Y “ C YpP + ^Yg^r 
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An. expression for (3 in terms of tbs controlled-airplane response and 
Pg is ' 


0 — Po "** Pg 

- - + (ja) 


The "block diagram of figure l4-(a) illustrates how the side-force and 
yawing -velocity feedbacks are included in dete rmining i|r/p g of the con- 
trolled airplane. 

The variable \jr is used as an input to the aileron servo in the 
type 2 autopilot. Therefore, an expression containing this autopilot 
coupling and the aerodynamic coupling between 0 and p is necessary 
to determine the 0/p g response. The aerodynamic coupling is included 
as was done for the calculation of the type 1 autopilot, with p deter- 
mined by the equation (Al). Figure 14 illustrates the manner of com- 
bining the transfer functions to determine the 0/p g response. 

In the determination of the roll response to rolling gusts 0/D0 g , 

the autopilot coupling between ty and 0 should be included as was done 
before in determining 0/p g . However, because of the very small aerody- 
namic coupling between 0 and i|f, it was believed that the t/D0 g 

response would be very small and that the yaw signal to the ailerons 
could be neglected. The response of 0 to D0 g is obtained from a 
single -degree-of -freedom calculation illustrated by the block diagram 
of figure 14(b). With a Value for 0/D0 g determined, it is possible 
to obtain i|r/D 0 g l 11 a manner similar to that used for the type 1 
autopilot. 


Type 3 Autopilot 

For the airplane in combination with the type 3 autopilot, the 
block-diagram "method of calculation for t/Pg is identical with that for 

the airplane in combination with the type 1 autopilot. Since the auto- 
pilot coupling between side force and aileron deflection is very strong 
in this case, it is recognized that the resulting large angle of roll 0 
might influence the motion in yaw response \|r at low frequencies. Wo 
method for including the roll effect could be devised; and, as a result, 
it is to be expected that these calculations will show more error in the 
lower range of frequencies than the block-diagram calculations used in 
the other cases. 
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In order to dete rmin e 0/p g , side force is determined as was done in 

the calculations for the type 2 autopilot, and this side force is used as 
an input to the aileron servo. The response of 0 of the uncontrolled 
airplane is determined from a consideration of a single degree of freedom 
in roll as was done in the case of the type 1 autopilot. The sum of the 
airplane response of 0 to Pg and the airplane response to the aileron 
deflection which resulted from the side-force input and the roll- 
stabilization feedback of the autopilot gives the controlled-airplane 
response. In determining the side force, 3 is assumed to be composed 
of -i|r, Pg, and the increment of 3 resulting from the lateral velocity 

of the airplane tha t arises because of the roll angle . The block diagram 
representing the analysis is shown in figure 15 (a). 

The simplified block-diagram method of calculating the responses to 
rolling gusts (fig. 15 (b)) is the same as for the airplane in combination 
with the type 1 autopilot. 
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TABLE I.- CHARACTERISTICS OF AIRPLANE USED IN CALCULATIONS 


Physical characteristics : 

Weight, lb 12,600 

Radius of gyration about X-axis kjj, ft 4.30 

Radius of gyration about Z-axis kg, ft 7*91 

Wing area, sq. ft 250 

Vertical tail area, sq. ft 55 

Wing span, ft 35*25 

Tail length, ft 15 

Distance from thrust line to center 

of area of vertical tail, ft 5*32 

Flight conditions: 

Airspeed, ft/sec 695 

Altitude, ft ■ 30,000 

Mach number 0.7 


Tyn nrHmenpH n nal data: 

H 

*x 

*Z 

% 2 ' 

Kg 2 


K xz 

c L 


50 

0.42 

0.112 

0.01485 

0.0504 

0.00062 

0.242 


Aerodynamic parameters: 


Side force Yawing Rolling 


CY Pwb 

. . -0.05 

Swb ’ ’ 

. . -0.055 

Cz 3 wb * * 

. . -0.11 

%t 

. . -0.53 

%. •• 

. . 0.175 

%t 

. . -0.05 

Cy • * 

x Dp 

. . 0 

Cj, . . 

“dp 

. . 0 

c 7 . . 

‘Dp 

. . 0 

%,* • • 

. . 0 

Cn • • 

^wb 

. . - 0.039 

C 7 * • 

P wb 

. . - 0.45 

° Y d 20 • * 

. . -0.045 

Cd D 20 * * 

. . 0.0147 

C 7 J . . 
Z D 2 < f > 

. . - 0.005 

Cy • • 

. . -0.050 

c n - * 

. . 0.0204 


. . -0.005 

*t 

^t 



° Y ^r • * 

. . 0.228 

C *$ T * * 

. . -0.077 

V ■ • 

. . 0.026 

C X • • 

. . 0 

^TLp. • • 

6 a 

. . 0 

%. ' ■ 

. . -0.086 

Cy • • 

i r 

. . 0.45 

Cn r 

. . -0.147 

C z . . 

r vb 

. . 0.004 





C 7 . . 

ir t 

. . 0.050 


Sidewash 

£|£ 0 



EASE H.- AIMOPXLOT GAH3, EQUIVAI£HT 8TABILITX 
DHOTAIIVES, AMD FOH4ULAS 


Autopilot galna 

Kqulralant stability darlratlToa 

Far all throe autopilots 

Budder servo, K r - 0.0322 radiana/volt 

Harm 

Aileron aervo, ■=> 0.0555 radlana/volt 

Mono 

Ky => 520 volt/ radian/ spaa 

* 2C n^K^Kr “ -2*60 

“ 1^.5 volta/radlaa 

c i$ - - ^■ o6e 

lor type 1 autopilot 

Kf - J12 volts /radian 

- Cjj^K^K,. “ -0.77 

For typo 2 autopilot 

Ky m 2J TOlts/radlan . 

0 2^, " ■ — 0 . XI 

J£| - 11 V volts /radian 

% - » -0.5^1 

^1 = 109 Tolta/aondlTiannlooal force 

S ■ ■ 43,15 

For typo 5 autopilot 

Kp ■=■ l,l80 volta/nondijnwnsioral foroo 

% - “ - 3,2? 

■ 312 volta/radlan 

C H “ C l - 12.8 
C *I* ' WW - ^ 

0^ 0 “ -0*77 

Far typo 3 autopilot vhen heading ocsmanA 

la ueod to determine n/f 0 


% 0 “ c Ye r ( Ic ^ + “fW ’ 2,26 + 5,0OD 

“ C ne r ( K (f + nc t) Ic r “ -0.77 + 1.3S0D 


C U “ c l. K P C ^r_( K + + “tyK “ 12,8 + fil>2D 

mj °b °r 
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Lift 




Figure 3.- Sample block diagram illustrating method of calculation of 
gust response of the airplane-autopilot combination. 
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Figure 5 .- Response to gusts of the airplane with controls fixed 
M = 0.7 i altitude, 30,000 feet. 
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(a) Side gusts. 


(b) Rolling gusts. 


Figure 6.- Response to gusts of the airplane in combination with the 
type 1 autopilot. M = 0.7j altitude, 30,000 feet. 
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to, radians / span 

(a) Side gusts. (t>) Rolling, gusts. 


Figust 8.- Response to gusts of the airplane in combination with the 
type 2 autopilot. = ll^j M = 0.7j altitude* 30*000 feet. 
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(a) Side gusts. 


(b) Rolling gusts. 


Figure 9.- Response to gusts of the airplane in combination with the 
type 3 autopilot. M = 0.7j altitude * 30^000 feet. 
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('•h'i Atmlane in combination with (o) Airplane in combination with 
"type 1 autopilot. type 2 autopilot. = 2J. 

Figure 10.- Continued. 
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(d) Airplane In combination vlth , (e) Airplane in combination with, 

type 2 autopilot. K^. = 114-. type 5 autopilot. 


Figure 10,- Concluded. 





.ding command for the airplane 
autopilots . . M =■ 0 . 7 j altitude 







Amplitude ratio 



Figure 12.- Variation of course with heading command for the airplane 
in combination with the three autopilots. M => 0.7j altitude, 
30,000 feet. 
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Airplane, free in roll only 



(a) Side gusts . 
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Airplane, free in roll only 
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(h) Rolling gusts. 


Figure 13.- Simplified transfer-function method of calculating gust 
response of airplane in combination with type 1 autopilot. 
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(h) Bolling gusts. 

Figure l4.- Simplified transfer-function method of calculat ing gust 
response of airplane in combination with type 2 autopilot. 
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